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ABSTRACT 

We performed a statistical analysis of half-power pulse-widths of the core compo- 
nents in average pulsar profiles. We confirmed an existence of the lower bound of the 
distribution of half-power pulse-width versus the pulsar period W50 ~ 2°. 45 P^"-^ 
found by Rankin (1990). Using our much larger database we found W50 = (2°. 51 ± 
0°.08) p-o.50±o.02 £qj, 21 pulsars with double-pole interpulses for which measurement 
of the core component width was possible. On the other hand, all single-pole inter- 
pulse cases lie in the swarm of pulsars above the boundary line. Using the Monte 
Carlo simulations based on exact geometrical calculations we found that the Rankin's 
method of estimation of the inclination angle a « asin{2°A5 P~'^'^ /W^q) in pulsars 
with core components is quite good an approximation, except for very small angles a 
in almost aligned rotators. 

Key words: stars: pulsars: general - stars: neutron - stars: rotation 



1 INTRODUCTION 

It is generally accepted that the elements of mean pu lsar pro- 
files c an be divided into core an d con al components (Rankinl 
l|l983l ). Rankin (1990, hereafter lR9oh . iRankinI j 19931 )). The 
core component usually occupies central parts of the profile 
and is often flanked by one or two pairs of conal compo- 
nents. The importance of the conal profile components was 
recently studie d statisti cally by Maciesiak, Gil & Ribeiro 
(2011; hereafter IPaper j ). These authors used recently large 
data sets of pulsar periods, profile widths and inclination 
angles resulted from large pulsar surveys conducted in re- 
cent years. They also compiled the largest ever database of 
pulsars with interpulses (IP), divided into the double-pole 
(DP-IP) and single-pole (SP-IP) cases (see Table [T] in this 
paper). The IP pulsars constitute about 3 percent (2 and 
1 percent for DP-IP and SP-IP cases, respectively) of the 
pulsar population. The DP-IP pulsars (almost orthogonal 
rotators) tend to be younger or middle aged, w hile SP-IP 
pulsars tend to be older (see Figure 6 in lPaper j) . This sug- 
gest a secular alignment of the magnetic axis towards the 
spin axis over a timescale of about 10^ years. 

The interpulse emission is important also for the statis- 
tics of the core profile components in pulsars. Rankin (1990; 
hereafter R90) measured the half-power (50% of the maxi- 
mum intensity) widths of core components in a number of 
pulsars (about 100), including those showing the interpulse 
emission. She found that the distribution of pulse-widths 
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has a clear lower bound 2°. 45 P""-^ populated by DP-IP 
cases, in which the inclination angle a must be close to 90°. 
The rest of pulsars including SP-IP ca ses li ed in a swarm 
of data points above the lower bound. IR9CI described this 
distribution by a simple mathematical expression 



Ws^ 



2°. 45 P~°-^/sina, 



(1) 



indicating that the core component widths depend only 
upon the pulsar period P and the inclination angle a be- 
tween the rotation and magnetic axes of the neutron star 
and they are almost insensitive to the impact angle /3 of the 
closest approach of the observer's line-of-sight to the mag- 
netic axis (for description of pulsar geometry see Paper I or 
Figure [Bll in the on-line Appendix [B]) . The above equation 
represents a useful method of estimation of the inclination 
angle a in pulsars with core component in the mean profile. 
In this paper we attempt to revisit and verify these impor- 
tant and intriguing r esults, u sing a more abundant database 
that we compiled in lPaper a) . 



2 IMPORTANCE OF THE CORE PULSAR 
EMISSION 

Table [T] contains 44 known IP cases divided into double- 
pole (DP-IP) and single-pole (SP-IP) cases (see IPaper j 
for details). This table presents for each case the values of 
VK50 predicted by Equation ((T| wi th a = 90° derived by 
IR90| . as well as those measured by IR90| and meas urements 
made in this paper based on IP studies presented in lPaper j . 
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[s] 2°.45P-0-5 


R90^ 


this paper 


DP-IP 




1. 


J0534 + 2200 


B0531 + 21 


0.033 


13.49 


13.5 ±1 




DP 


T/K 


2. 


J0627 + 0706 




0.476 


3.55 






DP 


KIO 


3. 


J0826 + 2637 


B0823 + 26 


0.531 


3.36 


3.38 ±0.1 


3.47 ±1 


DP 


T/K 


4. 


J0834 - 4159 




0.121 


7.04 




7.0 ±1 


DP 


K03 


5. 


J0842 - 4851 


B0840 - 48 


0.644 


3.05 




3.1 ±0.5 


DP 


Mil 


6. 


J0905 - 5127 




0.346 


4.17 




5.0 ±1 


DP 


WJ 


7. 


J0908 - 4913 


B0906 - 49 


0.107 


7.49 


7.5 ±0.4 


7.5 ±1 


DP 


T/K 


8. 


J1057 - 5226 


B1055 - 52 


0.197 


5.52 


<6 


5.8 ±1 


DP 


T/K 


9. 


J1126 - 6054 


B1124 - 60 


0.203 


5.44 






DP 


WJ 


10. 


J1244 - 6531 




1.547 


1.97 






DP 


K09 


11. 


J1413 - 6307 


B1409 - 62 


0.395 


3.90 




4.1 ±1 


DP 


Mil 


12. 


J1549 - 4848 




0.288 


4.57 




5.0 ±1 


DP 


A 


13. 


J1611 - 5209 


B1607 - 52 


0.182 


5.74 




5.1 ±1 


DP 


WJ 


14. 


J1613 - 5234 




0.655 


3.03 






DP 


R-MOl 


15. 


J1627 - 4706 




0.141 


6.52 






DP 


R-L 


16. 


J1637 - 4553 


B1634 - 45 


0.119 


7.10 






DP 


WJ 


17. 


J1705 - 1906 


B1702 - 19 


0.299 


4.48 


4.5 ±0.3 


5.5 ±1 


DP 


T/K 


18. 


J1713 - 3844 




1.600 


1.94 




2.1 ±0.5 


DP 


K03 


19. 


J1722 - 3712 


B1719-37 


0.236 


5.04 




6.0 ±1 


DP 


T/K 


20. 


J1737 - 3555 


B1734 - 35 


0.398 


3.88 




4.65 ± 1 


DP 


Mil 


21. 


J1739 - 2903 


B1736 - 29 


0.323 


4.31 




4.8 ±1 


DP 


T/K 


22. 


J1825 - 0935 


B1822 - 09 


0.769 


2.79 


2.8 


2.8 ±0.25 


DP 


T/K 


23. 


J1828 - 1101 




0.072 


9.13 






DP 


M02 


24. 


J1842 + 0358 




0.233 


5.08 




7.0 ±1 


DP 


R-L 


25. 


J1843 - 0702 




0.192 


5.59 




5.8 ±1 


DP 


H 


26. 


J1849 + 0409 




0.761 


2.81 




3.0 ±0.6 


DP 


L 


27. 


J1913 + 0832 




0.134 


6.69 






DP 


M02 


28. 


J1915 + 1410 




0.297 


4.49 






DP 


R-L 


29. 


J1932 + 1059 


B1929 + 10 


0.227 


5.14 


5.15 ±0.1 


5.3 ±0.2 


DP 


T/K 


30. 


J2047 + 5029 




0.446 


3.67 




3.88 ±0.5 


DP 


J 


31. 


J2032 + 4127 




0.143 


6.47 






DP 


C 


1. 


J0828 - 3417 


B0826 - 34 


1.849 


1.80 


40 


70 


SP 


T/K 


2. 


J0831 - 4406 




0.312 


4.39 






SP 


R-K03 


3. 


J0953 + 0755 


B0950 + 08 


0.253 


4.87 




30 


SP 


T/K 


4. 


J1107-5907 




0.253 


4.87 




50 


SP 


R-L 


5. 


J1302 - 6350 


B1259 - 63 


0.048 


11.21 




39 


SP 


T/K 


6. 


J1424 - 6438 




1.024 


2.42 






SP 


R-K03 


7. 


J1637 - 4450 




0.253 


4.87 




19 


SP 


R-L 


8. 


J1806 - 1920 




0.880 


2.61 




20 


SP 


M02 


9. 


J1808 - 1726 




0.241 


4.99 




97 


SP 


L 


10. 


J1851 + 0418 


B1848 + 04 


0.285 


4.59 




20 


SP 


T/K 


11. 


J1852-0118 




0.452 


3.64 




20 


SP 


H 


12. 


J1903 + 0925 




0.357 


4.10 




155 


SP 


Mil 


13. 


J1946 + 1805 


B1944 + 17 


0.441 


3.69 




7 


SP 


T/K 



'^ R90 lists 6 cases of DP-IP in he r Tab le 3 that she used to find a formal fit, but we were not able to confirm the 
interpulse in B033-45 (Vela pulsar) . |R90| in Table 2 lists few more DP-IP cases for which precise measurements of W50 
were not possible, so they were not included in her fit. Two additional cases B1055 — 52 and B1822 — 09 were taken 
from Tabic 5 in R90. |R90| listed only one case of SP-IP. 
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Among 31 DP-IP and 13 SP-IP cases 21 and 11 measure- 
ments of W50, res pectively, were possible. We utilized the 
ATNfQ data base (JManchester et alJlJOOsI ). where the val- 
ues of overall half-power widths were tabulated. If needed, 
the raw data for individual profiles were available for inspec- 
tion and/or processing]. Figure [T] shows the plot Wbo versus 
P for DP-I P cas e only. Red open symbols represent mea- 
surements of 'R9u presented in her Table 3, while black filled 
symbols represent our measurement^ presented in Table [T] 
in this paper. We added the case of longest period radio 
pulsar J2144-3933, as its W50 value (ATNF database) falls 
close to the value predicted by IR90| for DP-IP pulsars. Al- 
though, the interpulse has not been found in this pulsar, this 
seems consistent with its very narrow pulse profile (inter- 
pulse could be just missed by our lin e-of-sight). This p ulsar 
was identified as the core profile case (| Young et al.l 19991 ) and 
it had the narrowest pulse profile among all pulsars known 
at that time. 

The formal fit to all data points (excluding J2144 — 
3933) is 



W5o = (2°. 51 ±0°. 08) P' 



(2) 



which is represented by the black line in Figure [T] This 
is very close to the fit corresponding to the data of R90 
(marked in red) W50 = (2°.46 ± 0°.09) p-0-50±o.02 ^ ^^^^^_ 
sented by red line in Figure [T] Therefore, the lower bound 
found by IR9C|| still appears clearly in our larger data set 
analysed in this paper. This is illustrated in Figure (2] which 
contains W50 measurements of DP-IP (red), S P-IP (green) 
and core components of other pulsars (data of |R90| marked 
by black dots (Table ICll in the on-line Appendix [C]) and 
our new measurements marked by black stars (Table [C2l in 
the on-line Appendix [C|. At least from the first sight the 
Rankin's Equation ([l]) seems well established. The impor- 
tant case of extremely broad profile pulsar B0826 — 34, with 
core component width being close to 70°, is marked (see also 
caption to Figure [5]) . 

Figure [2] was prepared using W50 measurements of the 
core profile components. However, one can ask the question 
how would this figure look like if one used all W50 measure- 
ments available in the ATNF database (without distinctions 
into core and conal components). Some of the measurements 
corresponding to profiles dominated by prominent core com- 
ponents would just reproduce the Figure [51 However, those 



^ http://www.atnf.cisiro.au/rcsearch/pulsar/psrcat/ 
^ http://datanet.csiro.au/dap/ 

^ We measured or estimated VF50 of the core components in 
21 DP-IP and 11 SP-IP pulsars. In a few cases of DP-IP cases 
it was possible to obtain the so-called log-polar representation 
IIHankins &: Fowled [1986) of the profile and measure the core 
widths quite accurately. However, in most cases the estimate was 
done by means of the visual inspection of the profile and a rea- 
sonable judgement. For example, if the profile was apparently 
complex but its components could not be resolved we used the 
following method: we estimated the half-power width of the en- 
tire envelope of main pulse or inter-pulse and then we divided the 
obtained result by an integer following from the complexity of the 
envelope (mostly 3). In such cases quite conservative e rrors were 
given. Our measurements were consistent with those of lRQCl in all 
cases where the comparison was possible. In fitting procedures we 
used mostly our new measurements except five cases for wich the 
estimates given by R90 were more accurate than ours. 



)f„=(2°.51±0.08) p-0-50±0.02 




Period P [s] 



Figure 1. Plot of W50 vs. P for 21 DP-IP pulsars. Five cases 
marked by red circles were taken from Table 3 in R90 and the rest 
of points represent our new measurem ents ( see Table [TJ . The red 
solid line represents the original fit of|R90|, while the black line 
represents the fit to all data points. We also added the longest 
period radio pulsar J2144 — 3933 marked in blue to show that it 
follows the boundary. However it was not included in the fitting 
procedure as it did not had the IP emission. 



100 



,B0826-34 




Period P [s] 



Figure 2. Similar as in Figure [T] but with measurements of 
core widths for a number of other pulsars, including S P— IP cases 
marked in green. Dots represent the measurements of [R90| while 
stars represent our measurements. For DP-IP cases mar ked in 
red, the open symbols represent measurements of IR9C|| . while 
filled symbols represent our measurements. The important case 
of B0826 - 34 is also marked. 
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Figure 3. Plot of pulse-widths W^o measured at 1.37 GHz for 1522 normal pulsars (excluding millisecond and other recycled pulsars; 
see Paper I for details) from ATNF database, with added 21 cases of DP-IP (marked in red) and 11 cases of SP-IP (marked in green) 
taken from Figure [2] (Table [TJ- The black solid line represents the lower limit from Figure [2] rescaled to the frequency 1.37 GHz. The two 
dashed grey lines are marked to illustrate the scatter of red data points and their errors around the black solid line W50 = 2°. 37 P""'^. 
The case of J2144 — 3933 is marked by blue cross, while blue circles correspond to pulsars with partial cone profiles. The orange stars 
represent RRATs, the magenta stars represent young energetic pulsars and green stars represent pulsars discovered as strong single pulses 
(see Table [2] for the list of special cases) . 



measurements which correspond to the conal profiles should 
lie above the Lower Boundary Line 2°. 5 P~'^-^ (LBL here- 
after), which should still appear as the prominent feature. 
Figure [3] presents a plot of all 1522 values of W50 versus 
pulsar period, which were available in the ATNF database 
(excluding millisecond and other recycled pulsars). Appar- 
ently, the bulk of the measurements lies above the LBL 
(W50 = 2°.37P-°'^ after rescaling to 1.37 GHz). The widths 
of DP-IP (red symbols with error bars) and SP-IP (green 



symbols) were added for the reference and comparison with 
Figure [1 

It is natural why the absolute majority of points lies 
above the lower boundary line. Indeed, if these are core 
components and/or conal components surrounding the core 
ones, then they must obey the lower boundary limit as dis- 
cussed in this p aper (for geo metry of conal compo nents see 
iRankinI (|l993l ). IGKS93I and iKramer et aP (| 19941 )). There- 
fore, the real question is why there is a small number of 
points below the LBL (assuming that these measurements 
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Table 2. Special pulsars appearing in Figure [3] near or below the 
lower boundary linef. 



Name J 


P[s] 


H/«S0[dcg] 
1.4GHz 


W5o[deg] 
1.4GHz 


Ref. 






Partial cone 




J1759 - 2205 


0.461 


3.61 


3.05 


MRU 


J1912 + 2104 


2.232 


1.64 


1.60 


MRU 


J1917 + 1353 


0.194 


5.56 


7.40 


MRU 


J1941 - 2602 


0.403 


3.86 


3.66 


MRU 


RRATs 


J0735 - 62 


4.865 


1.11 


0.15 


KKLll 


J1047 - 58 


1.231 


2.21 


1.08 


KKLll 


J1226 - 32 


6.193 


0.98 


0.70 


KKLll 


J1423 - 56 


1.427 


2.05 


1.14 


KKLll 


J1514 - 59 


1.046 


2.40 


1.14 


KKLll 


J1554 - 52 


0.125 


6.93 


2.88 


KKLll 


J1654 - 23 


0.545 


3.32 


0.66 


KKLll 


J1707 - 44 


5.764 


1.02 


0.76 


KKLll 


J1724 - 35 


1.422 


2.05 


1.49 


KKLll 


J1807 - 25 


2.764 


1.47 


0.52 


KKLll 


J1819 - 1458 


4.263 


1.19 


0.25 


HEYll, KHS09 


J1826 - 1419 


0.771 


2.79 


0.93 


MLK09 


J1841 - 14 


6.598 


0.95 


0.14 


KKLll 


J1854 + 03 


4.558 


1.15 


1.25 


KKLll 


J1848 - 12 


6.795 


0.94 


0.11 


MLL06 


J1913 + 1330 


0.923 


2.55 


0.78 


MLK09 






Young energetic 


PSRs 




J0531 + 2200 


0.033 


13.49 


4.4 


R90 


J1028 - 5819 


0.091 


8.12 


0.63 


KJK08 


J1225 - 6035 


0.626 


3.10 


2.24 


KBM03 


J1855 + 0527 


1.393 


2.08 


0.78 


KEL09 


Single Pulse 


J0627 + 16 


2.180 


1.66 


0.05 


DCM09 


J1909 + 06 


0.741 


2.85 


0.73 


DCM09 


J1946 + 24 


4.729 


1.13 


0.31 


DCM09 



t There were also two recycled pulsars (J1410 — 7404 and J2129 + 
1210G) below this line, which we did not show according to strat- 
egy adopted in Paper I. 

are not erroneous)? Taking into account the scatter of red 
data points and their errors in Figure[3]it is instructive to in- 
troduce the Lower Boundary Belt W50 = (2°.37+Jo;g?)P"°-^ 
(LBB hereafter; marked by two grey dashed lines), which 
should not be mistaken with the formal fit expressed by 
Equation ((2]). We took every effort to identify points ly- 
ing below or within the LBB (Table [2]). ft is extremely in- 
teresting that most of these points were Ro tating RAdio 
Transients (RRATs; iMcLaughhn et all (|2006l )') marked by 
orange stars, or other sporadically emitting pulsars. RRATs 
are repeating sources with pulsar underlying periodicity, de- 
tectable via their extremely strong sporadic single pulses. A 
good example is the case of Jf8f9 — f458 marked by the 
orange error bar in Figure [31 whose length represents in- 
dividual measurements of different events (Hu et al. 20f f), 
with the average value of W50 = 0°.25 (Table ^. Inter- 
estingly, the width of the mean profile obtained by folding 
of large number of individual events of this relatively fre- 
quent RRAT is about fO°, well above the LBB. However, 



the most interesting fact is that the width of the central 
co mponent of this profile is about f°.4 ± 0°.4 (Figure L 
in iKarastergiou et al.l (|2009l ')'). which fits to the LBL/LBB 
W50 = 2°. 37 (4.3)-°'^ = 1°.14 very well. This strongly sug- 
gest that most of data points below the LBB in Figure[3]rep- 
resent a single events of RRATs activity, whose pulse-widths 
are much smaller than their underlying pulse windows. 

Another small group of three data points below the LBB 
represents the Single Pulse (SP) det ected pulsars (gree n 
stars) in the Arecibo PALFA survey (JDeneva et al.l 120091 '). 
One of them J0627 -I- f 6 has the narrowest ever measured 
width of 0°.05 (Table [2]), which corresponds to the lowest 
known duty cycle of about 0.01 percent (compared to the 
typical value of about 5 percent). 

Yet another group of data points below the LBB rep- 
resents the young energetic PSRs (Table (2]) marked in ma- 
genta. It is interesting that the main pulse of the Crab pulsar 
(marked by diamond) seems to belong to this category as 
well. It is well known that the precursor of the main pulse 
in the profile of B0531 + 21 behaves like typical core compo- 
nent, while the main- pulse and the interpulse show rather 
untypical properties (e.g. R90|). 

We also checked the group of pulsars with partial cone 
profiles but no single case was found below the LBB in Fig- 
ure[3l In the long list o f such pulsars recently published by 
iMitra &: RankinI l|201ll 'l we identified only four cases (Table 
[2} lying within the LBB (blue circles). This paradoxical at 
first sight result means that the opening angles of core and 
conal beams are about the same. 

In summary, the lower boundary with the slope follow- 
ing p-^^-'' is a prominent feature in the distribution of pulse- 
width 14^50 versus pulsar period P, especially if extended to 
the concept of narrow lower boundary belt (Figure ^ with 
the width determined by scatter and errors of W50 mea- 
surements of core components in DP-IP pulsars (red data 
points). There is a small number of data points below the 
LBB but all of them represent some unusual species, mostly 
transients and intermittent pulsars, whose sporadic emission 
is much narrower than their underlying pulse windows. Fi- 
nally, it is worth noting that the existence of LBL or even 
LBB is not depending neither on the method^ nor on the 
accuracy of the measurements of core widths in DP-IP pul- 
sars. Indeed, the LBL is still apparent when the red symbols 
are removed from Figure [31 as clearly visible in Figure [5] 



3 GEOMETRY OF THE CORE EMISSION 

The period dependent lower boundary in W50 pulse-widths 
revealed in pulsar data, translates into the period dependent 
lower limit for the half-power beam widths of the transversal 
structures in the overall pulsar beam. Geometry of pulsar 
emission is described by the well known formula 



W50 = 4 arcsin 



sin [{p5o + P) /2] sin [(pso - P) /2] 



11/2 



sin a sin{a + 0) 



, , (3) 

(JGiil l|l981^ : lGil■ Gronkowski fc Rudnickil (|l984h . Gil, Kijak 
& Seiradakis (1993); hereafter GKS93), where a is the in- 
clination angle (of the magnetic axis to the rotation axis), 
P is the impact angle (of the closest approach of the line- 
of-sight to the magnetic axis) and pso is the opening angle 
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of the pulsar beam corresponding to half-power maximum 
beam width (FWHM). This equation holds generally on the 
sphere centered on the pulsar and the only assumption about 
the pulsar beam that has to be satisfied is the phase symme- 
try about the fiducial plane fo containing both pulsar axes 
and the observer's direction (Figure IBll in the on-line Ap- 
pendix B). One should realise that the illustration of circular 
beam presented in Figure [Bll represents only a cross-section 
of a certain intensity level of the beam with the pulsar cen- 
tered celestial hemisphere. It says nothing about the internal 
structure of the beam, neither the angular nor the radial in- 
tensity distribution. The core pulsar beam can be modelled 
by a gaussian distribution of the mean intensity 7, symmet- 
ric about the magnetic axis (beam axis). Thus, it can be 
described by the function 



7 = 7o exp[-fc(p/po)^], 



(4) 



where 7o = 7(/3 = 0) = Imax, I{po) = 7o exp(-fc) = Imin 
and k is an arbitrary steepness parameter. For a certain 
impact angle « /3 < p the observer will detect a local 
intensity maximum 

7™(/3) = 7o exp[-fc(/?/po)'] 
and thus we can write Equation Q in the form 



7(p,/?)=7„(/J)exp 



^ip'-P') 



Po 



(5) 



where Im{P) is the local peak intensity which can be de- 
termined from observations. R90 used slightly different ap- 
proach and assumed that the angular distribution of the 
core pulsar beam is well approximated by a bivariate gaus- 
sian function 

7(/?, ^) ex exp - [(/? - A))V2P0 + (^ - Vof/2p%], 

whose breadth (FWHM) is invariant along any set of par- 
allel paths. Thus, she concludes, "the angular width of core 
beams can only be weakly dependent on the centrality of 
traverse and the effect of the impact angle /? on the ob- 
served width of the core components may be quite weak" . 
However, one should realise that approximation described 
by the above equation is a true representation of gaussian 
intensity distribution of the pulsar beam only in a special 
case when tp r^ p, that is /3 f» 0° and a ~ 90° (Equation 
[Sjl. Only then the above function I{/3,ip) is equivalent to 
I{p,l3) expressed by our exact Equation ([Sj. In general, al- 
though the projection of the line-of-sight cut through the 
gaussian beam is indeed independent of the impact angle /3, 
the observed pulse-width (not a beam width) W = 2ip may 
depend on both a and /3. To examine this problem exactly 
we will use Equation ([3]) equipped with the relationship 



P50(/3) 



I 2 , q2\1/2 



(6) 



valid only for a gaussian beam shape (see Figure [B2l for il- 
lustration). Here pso = P5o(/3 = 0) is the half-power beam 
width, while psoiP) is the opening angle of pulsar beam mea- 
sured at FWHM level by the observer with the impact angle 
p. Following Figures [T] and [2] it is quite reasonable to adopt 
for the half-power width of the core beam 



P50 



^o 2p-l/2^ 



(7) 



which seems quite well established observationallj|j. The 
natural P~^" factor undoubtedly means that the radiation 
associated with W50 of core components is emitted tangen- 
tially to dipolar magnetic field lines, but the value of the 
numerical factor close to 1°.2 is not yet clear (see Conclu- 
sion section). 

The results of calculations of W5o(a, /3, p(/?)) for differ- 
ent values of the pulsar period P are presented in Figure 
13] Three different panels correspond to different periods P 
equal to 1, 0.1 and 0.01 seconds, respectively. Four solid lines 
in each panel correspond to different values of the impact /3 
and opening p5o(/3) angles marked by red, green, dark blue 
and light blue colours (the latter corresponds to the value 
of about 0.9 p5o(/3)). The dashed orange line represents the 
maximum difference between the solid lines (red and light 
blue). As one can see for P = Is a significant residual occur 
only for very small inclination angles a < 5° (almost aligned 
rotators). For shorter periods the residuals are larger. How- 
ever, even for 10 millisecond periods the errors made due to 
neglecting /? are negligible for inclinations angles a > 30°. 
Thus, one can generally state that neglecting /3 in estimat- 
ing the core width is quite well justified by results of strict 
calculations, especially for larger inclinations angles. For or- 
thogonal rotators (DP-IP) it can be treated as an exact 
result. Thus, Equation ([l]) proposed bv IR9C1 should be re- 
garded as a very good tool for estimating the inclination 
angles using the core pulse-widths. In Figure \S\ we super- 
posed on the data from Figure [2] the results of simulations 
of W50 performed using a technique developed in Paper I ap- 
plied to core emission (see on-line Appendix 1^. Although 
these were strict calculations, including simulations of both 
inclination angle a (different ranges of this angle are marked 
with different colours) and impact angle P, their results are 
fully co nsiste nt with Equation ([ij. This equation was de- 
rived bv IR90| assuming /3 = 0, based on the claim that the 
core width is almost completely insensitive of the value of 
the impact angle /3. Pulsars with interpulses corresponding 
to DP-IP and SP-IP models are marked by red and green 
symbols, respectively. All DP-IP cases lie on the boundary 
line (to within error bars), as expected. A very interesting 
case of almost aligned rotator PSR B0826 — 34 is marked in 
Figure [S] The inclination angle of this extremely broad pro- 
file pulsar was estim ated as having v ery small value a < 5 ° 
(iGuptaet al.ll2004l ) or even a < 1° tesamdin et al.ll2005l ). 
which is consistent with results of our simulations. 



4 INCLINATION ANGLES AND 
PULSE-WIDTHS 

As demonstrated in Figure 4 using Monte Carlo simulations 
we can calculate the pulse-widths W50 of the core emission 
for large population of pulsars. These calculations are based 
on Equations ([iJ and (|3}, with a and /3 values drawn from 
suitable distributions described in Paper I (see on-line Ap- 
pendix A). On the other hand we can calculate values of W50 
using Rankin's Equation ([1} , with a and P values generated 

* The exact value that should appear in Equation (0 is just a 
half of the 2.45 factor from Equation JTJ. Of course such a four 
significant figures accuracy is not important, we nevertheless used 
it in our numerical calculations. 
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Figure 4. Calculations of W50 vs. a using Equations ((Sjl, lO and 
JT)!- Different panels correspond to different periods P, while four 
different colours of solid lines correspond to different values of 
j3 and p5o(/3). The asymptotic values 2p5o(0) are marked in the 
lower right corner of each panel. The dashed orange lines describe 
the maximum difference between the solid colour lines. 



by Monte Carlo simulations. Figure [B] presents comparison 
of W50 values obtained by these two different methods. Six 
different ranges of a values are marked by different colours, 
while values of /3 are represented by the length of vertical 
bars with the same scale as that of the ordinate axis. Of 
course, the correlation is perfect for very small /3 ~ 0°. As 
one can see the correlation is also very good for larger /3 
(corresponding to longer vertical bars), provided that in- 
clination angles are not too small (a > 20°). The points 
above and below the perfect correlation line correspond to 
positive and negative values of the impact angle /?, respec- 
tively. Generally one can say, that except of cases of almost 
aligned rotators both methods give very similar results. This 
further confirms usefulness of Rankin's method expressed by 
her simple Equation ((T}. 

Finally, we can make a direct comparison of Rankin's 
method of inclination angle estimations represented by 
Equation ([ij with strict calculations based on Equations (|3]) , 
([6| and (O ■ This comparison is presented in a form of scat- 
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Figure 5. The results of the Monte Carlo simulations of pulse- 
widths superposed on the Figure[2] Six different ranges of the in- 
clination angle a are marked by different colours. The DP-IP and 
SP-IP cases are marked by red and green symbols, respectively. 
The opened red symbols represent measurements of R90, while 
filled red symbols represent our new measurements (see Table 
[TJ. The solid line represents the fit to red marked data points ex- 
pressed by Equation ^. For a very wide profile pulsar B0826 — 34 
with SP— IP the simulations indicate the value of the inclination 
a ngle a < 1°, wh i ch is consistent with an independent estimate 
of lEsamdin et al.l l|2005l ) o = 0°.5 for this almost aligned rotator. 
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Figure 6. Comparison of W50 values obtained by means of Equa- 
tion JIJ and simulations based on Equations (|3]l, ^ and JT]!. Six 
different ranges of a values are marked by different colours, while 
values of /3 are represented by the length of vertical bars with the 
same scale as that of the right-hand side ordinate axis. 
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Figure 7. Scatter plot of inclination angles a obtained using 
Rankin's method (Equation JTJ) and Monte Carlo simulated val- 
ues. Red and green colours on the scatter plot correspond to pos- 
itive and negative values of the impact angle /3, respectively. The 
blue solid line represents the maximum relative difference between 
two methods (in percents according to the right-hand side ordi- 
nate scale). 



ter plot in Figure [T] The ordinate (vertical) axis represents 
a calculated from Equation ([ij with W50 calculated from 
Equations (|3}, ([6} and ([7}, where a, P and P were taken 
from Monte Carlo simulations, while the abscissa (horizon- 
tal) axis represents directly simulated values of a. Red and 
green colours correspond to positive and negative values of 
P (for details of simulations and detection conditions see 
Appendix IX} . The blue line represents the maximum rela- 
tive difference between scattered points (in the same scale 
as for a axes but expressed in percents), thus it describes an 
accuracy of determination of a. As one can see, for almost 
orthogonal rotators the difference between the two methods 
is about 10 percent. However, for almost aligned rotators the 
error can reach or even exceed 100 percent. For intermediate 
values of a the accuracy is in the range of 20 - 30 percent. 
This means that Rankin's method of estimating the incli- 
nation angle from the width of the core component is quite 
good an approximation, provided that these values are not 
too small. In particular, this method should not be applied 
to extremely broad profile pulsars for which one can expect 
very small inclination angles a ^ 5° . 



5 CONCLUSIONS 

The most important result of our paper is a clear confir- 
mation of the lower limit in the distribution of half-power 
pulse-width W50 in core components of pulsar profiles when 
presented versus the pulsar period P using much bigger 
database than that used originally by IR90| . As illustrated 
in Figures [T] [5J |3] and |S] this limit is represented by the 
line W50 = 2°.5P-°-^ at 1 GHz and W50 = 2° .37 p-°'^ 
at 1.4 GHz. It seems that more appropriate description of 
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Figure 8. Plot of pulse-width W50 versus period P for 1495 
normal pulsars from the ATNF database. 



the lower boundary limit is the lower boundary belt (LBB) 
W50 = (2°.37+oo;87) P'°-^, which accounts for scatter and 
errors in W50 measurements. Strictly speaking, the widths 
of the LBB was determined by the scatter and errors of 
DP-IP cases marked as red symbols in Figure 3. However, 
the normal data points seem to populate the belt marked 
in Figure 8 in similar manner as the DP-IP cases do in Fig- 
ure 3. The concept of LB B is b ased not only on core widths 
of DP-IP pulsars (like in IR90| ) but here it is based on all 
W50 measurements of normal pulsars available in the ATNF 
database, without distinctions into core and conal compo- 
nents, of course. This is dramatically illustrated in Figure [8] 
containing 1495 data points. This figure is a reproduction of 
Figure[3]with all special points described in the legend omit- 
ted. Apparently, the existence of the lower boundary limit in 
the pulse-width measurements discovered by Rankin (1990) 
does not depend neither on the method nor on the errors 
of measurements of the core widths in DP-IP pulsars. The 
unweighted formal fit to the points within LBB (between 
dashed lines) is l^so = (2°. 86 ± 0°.07) p-0-5i±002. 

The existence of the lower bound in the distribution of 
W50 versus pulsar period P is very intriguing and requires 
solid and convincing explanation. The exponent —0.5 is ev- 
idently related to the geometry of dipolar field lines in the 
radio emission region, while the value o f the numerical factor 
about 2°. 5 is a little more mysterious. lR9d argued that the 
core component widths are intimately related to the polar 
cap geometry at the stellar surface, whose bundle of the last 
open dipolar field lines have an angular extent (opening an- 
gle 2p) very close to 2°.45P~^'^. Although this natural and 
very appealing explanation is commonly accepted, it is not 
free from theoretical and interpretational problems. First of 
all, the coherent radio emission cannot originate at or very 
near the polar cap, because there are no plasma instabili- 
ties available in this region. Recently, a modern version of 
this idea is discussed, with arguments that the core emission 
originates below the conal one but the former is emitted well 
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above the polar cap surface. For example. iMitra et al.1 (|2007l ) 
argued that in PSR B0329 + 54 the emission of its core com- 
ponent originates over a range of alt±tjrdes_(up to about 100 
km) above the star surface (see also I Gill ([l99lj)). This broad- 
ens the required opening angle by a factor of about square 
root of the normalised emission altitude, that is about 3. If 
this is a case then one has to conclude that the core radio 
emission is related to an innermost subset of the open dipo- 
lar magnetic field lines. The problem therefore is to identify 
structures within the overall pulsar beam with an angular 
extent corresponding to the observed lower boundary limit. 
We plan to devote to this problem the forthcoming Paper 
III of this series. 

The interpretational problems mentioned above do not 
affect the approximate method of estimating the inclination 
angle a in any pulsar with a core component proposed by 
[R90. Wc have shown by means of Monte Carlo simulations 
based on the strict geometrical calculations that this method 
offers quite good an approximation with the accuracy of 5, 
10, 20 and 50 percent for q ~ 90°, a ~ 50°, a ~ 25° and 
a ~ 10°, respectively (Figure [7]). This method cannot be 
used in almost aligned cases with very small values of a. 
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ON-LINE MATERIALS 



APPENDIX A: MONTE CARLO SIMULATIONS 
OF THE CORE EMISSION 

To simulate the core pulsar emission we use the technique 
described in our Paper I and adapt it to the specifics of 
coral beams. We will simulate only normal pulsars, that is 
excluding all millisecond, binary and other recycled pulsars. 
In addition to Equation ^ we will use Equations © and 
(O valid exclusively for gaussian core beam. Following a 
technique used in Paper I we will usually make a number of 
subsequent steps. 

1. Generate the pulsar period P as a random number with 
the parent probability density function /(P) corresponding 
to log-normal distribution function (Equation (8) in Paper 
I, which seem to best suited statistically for representing the 
true parent distribution of pulsar periods). 

2. Generate the inclination angle a as a random num- 
ber with the parent probability density function /(a) corre- 
sponding to Equation (7) in Paper I, which seems to be best 
suited for reproducing the right ratios of interpulse cases in 
pulsar population. 

3. Generate the observer angle ^ = q -f /3 as a random 
number with the parent probability density function /(^) = 

4. Calculate the impact angle /? = ^ — a. 

5. Calculate the opening angle pso for the 50 percent of 
the maximum intensity of the core component according to 
Equation ((6]). 

6. Check the detection condition: pso > 3°P~^". To the 
next steps only detected pulsars are taken. The numerical 
factor 3° for low intensity edge of gaussian core beams was 
adopted arbitrarily, but its value influences only the num- 
ber of detected pulsars and has no important meaning for 
conclusions of this paper. 

7. Calculate the pulse-width of the core component at 
the 50 percent of the maximum intensity level: 

a) W§^^ according to the Equation ([T|) using simulated 
a and P values. 

b) W^5o™ according to the Equation ((3]), where pso is cal- 
culated from Equations © and ((7]). 

8. Calculate a^^" angle according to Equation ([l]), where 
VFso = Wit- 

For more technical details on Monte Carlo simulations 
of pulsar radio emission see Paper I. 




direction to 
observer 



Figure Bl. Geometry of the pulsar radiation. Schematic picture 
of the celestial hemisphere centered on the pulsar. The fiducial 
plane </> = is set by the rotation fJ and the magnetic m axes. The 
following angles are marked: longitudinal phase angle measured 
from the fiducial plane (^ = 0, the inclination angle a, the impact 
angle /?, the observer angle £^ = a + jS, the opening angle of the 
beam of radiation p and the polarisation angle i/i. 



APPENDIX C: NEW DATABASE OF 
PULSE-WIDTH Wtm OF CORE COMPONENTS 

In this section we present a number of tables containing 
measurements of FWHM pulse-widths of core components 
that were possible to make judging from published profiles. 
The references to appropriate papers are given. 



APPENDIX B: 
RADIATION 



GEOMETRY OF PULSAR 



Figure iBll represents a general view of the geometry of pul- 
sar radiation emitted tangentially to dipolar magnetic field 
lines and intersected with the celestial sphere centred on the 
neutron star. 

Figure IB2I presents 3-D geometry of the gaussian core 
pulsar beam, while Figure IB3I presents intersection of this 
gaussian beam with the pulsar centered celestial sphere. 
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Figure B2. The scheme of the gaussian-shaped distribution of 
the radiation intensity on the pulsar beam centred on the mag- 
netic axis m inchned to the rotation axis n at a angle. Two 
different trajectories are marked: central (red) and side (green), 
corresponding to the impact angle /? = and arbitrary 0, respec- 
tively. In both cases the 50 percent level of the local maximum 
is showed. It is easy to see that in both cases the half-width of 
the appropriate 2-D Gauss function is the same (pso). However, 
the corresponding widths of the observed radiation are different 
050 (/3) < </'5o(0). Note that p|q(/3) = P^q + P^ for gaussian beam. 




Figure B3. The gaussian-shaped pulsar beam presented in Fig- 
ure IB2I is projected on the spherical surface centered on the 
NS. Note that W5o(/3) < W50, where WsoijS) = 2<j>5o{l3) and 
W50 = 2050. 
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Table CI: Pulse-widths of the core components at the 50% maximum 
intensity level taken from Rankin (1990). 



Pulsars with interpulse 


No. 


Name B 


VK50 n 


P[s] 


No. 


Name B 


W,M ° 


P s 


1 


0531+21p 


13.5 


0.0331 


6 


1055-52 


5.3 


0.197 


2 


0531+21 


4.4 


0.0331 


7 


1702-19 


4.5 


0.2989 


3 


0823+26 


3.38 


0.5306 


8 


1822-09 


2.8 


0.7689 


4 


0826-34 


40 


1.8489 


9 


1929+10 


5.15 


0.2265 


5 


0906-49 


7.5 


0.1067 


10 


1736-29 


- 


0.3228 


Pulsars without interpulse 


Pulsars St type (Single) 


1 


0105+65 


5. 


1.2836 


24 


1844-04 


6. 


0.5977 


2 


0136+57 


7. 


0.2724 


25 


1859+03 


5.25 


0.6554 


3 


0154+61 


6. 


2.3517 


26 


1900+05 


6.8 


0.7465 


4 


0355+54 


8. 


0.1563 


27 


1900+01 


4.4 


0.7293 


5 


0540+23 


8.6 


0.2459 


28 


1907+02 


4.1 


0.9898 


6 


0611+22 


7.5 


0.3349 


29 


1907+10 


6.1 


0.2836 


7 


0626+24 


7.2 


0.4766 


30 


1907-03 


6.3 


0.5046 


8 


0740-28 


10. 


0.1667 


31 


1911-04 


3. 


0.8259 


9 


0835-41 


3.7 


0.7516 


32 


1913+10 


4.3 


0.4045 


10 


0940-55 


7. 


0.6643 


33 


1914+13 


6. 


0.2818 


11 


0942-13 


4.6 


0.5702 


34 


1915+13 


6. 


0.1946 


12 


0959-54 


3.9 


1.4365 


35 


1924+16 


5.7 


0.5798 


13 


1154-62 


13. 


0.4005 


36 


1929+20 


5.3 


0.2682 


14 


1240-64 


7.2 


0.3884 


37 


1930+22 


6.5 


0.1444 


15 


1449-64 


8.5 


0.1794 


38 


1933+16 


5.25 


0.3587 


16 


1556-44 


9.0 


0.2570 


39 


1946+35 


5.5 


0.7173 


17 


1557-50 


7.0 


0.1926 


40 


1953+29 


40 


0.0061 


18 


1641-45 


6.7 


0.4550 


41 


1953+50 


5.8 


0.5189 


19 


1642-03 


4.2 


0.3876 


42 


2002+31 


2.24 


2.1112 


20 


1706-16 


5.5 


0.6530 


43 


2053+36 


9.35 


0.2215 


21 


1749-28 


5. 


0.5625 


44 


2113+14 


7.4 


0.4401 


22 


1839+09 


4.5 


0.3813 


45 


2217+47 


5. 


0.5384 


23 


1842+14 


8.2 


0.3754 


46 


2255+58 


10. 


0.3682 






Pu 


Isars T ty 


pc(Tr 


iple) 






1 


0329+54 


5.5 


0.7145 


20 


1907+00 


2.43 


1.0169 


2 


0450+55 


7.2 


0.3407 


21 


1907+03 


15. 


2.3302 


3 


0450-18 


7. 


0.5489 


22 


1911+13 


4.4 


0.5214 


4 


0656+14 


8. 


0.3848 


23 


1913+16 


14. 


0.0590 


5 


0736-40 


14. 


0.3749 


24 


1914+09 


6. 


0.2702 


6 


0919+06 


5. 


0.4306 


25 


1916+14 


3.5 


1.1810 


7 


1112+50 


4.6 


1.6564 


26 


1917+00 


2.18 


1.2722 


8 


1221-63 


6. 


0.2164 


27 


1919+14 


8. 


0.6181 


9 


1451-68 


11.5 


0.2633 


28 


1920+21 


3.4 


1.0779 


10 


1508+55 


5. 


0.7396 


29 


1926+18 


3.9 


1.2204 


11 


1541+09 


32. 


0.7484 


30 


1952+29 


4.8 


0.4266 


12 


1558-50 


4. 


0.8642 


31 


2003-08 


10. 


0.5808 


13 


1604-00 


5.1 


0.4218 


32 


2020+28 


4. 


0.3434 


14 


1700-32 


4. 


1.2117 


33 


2045-16 


4.0 


1.9615 


15 


1727-47 


5. 


0.8298 


34 


2111+46 


11.7 


1.0146 


16 


1747-46 


3.5 


0.7423 


35 


2224+65 


9.5 


0.6825 


17 


1804-08 


7. 


0.1637 


36 


2319+60 


6. 


2.2564 


18 


1821+05 


5.4 


0.7529 


37 


2327-20 


2. 


1.6436 


19 


1826-17 


9. 


0.3071 
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Tab. 


Cllcont. 














No. 


Name B 


w,o n 


P[s] 


No. 


Name B 


wso n 


P[s] 


Pulsars M type (Multiple) 


1 


0523+11 


4.2 


0.3544 


8 


1831-04 


23. 


0.2901 


2 


0621-04 


6. 


1.0390 


9 


1857-26 


8. 


0.6122 


3 


1039-19 


4. 


1.3863 


10 


1905+39 


4. 


1.2357 


4 


1237-F25 


2.7 


1.3824 


11 


1910+20 


3. 


2.2329 


5 


1737-^13 


4.1 


0.8030 


12 


1919+21 


3. 


1.3373 


6 


1738-08 


4. 


2.0430 


13 


2028+22 


3.5 


0.6305 


7 


1745-12 


4. 


0.3941 


14 


2210+29 


4.3 


1.0045 



Table C2: New pulse-widths of the core components W50 obtained in this 
paper. The pulse-widths W at different frequencies / were corrected to 
about 1 GH z accordi ng: to W oc p oc r^' (/) oc /~ ' (for details see 
footnote 4 in lPaper J) . 



WEISBERG 1999 - 1.418 GHz 


Name B 


P 

[s] 


1^50 (1.418 GHz) 

n 


WrT^ (1 GHz) 

n 


0609+37 


0.298 


6.2 


6.4 


0940+16 


1.087 


9.1 


9.4 


1848+04 


0.285 


18.6 


19.3 


1859+01 


0.288 


4.6 


4.8 


1904+06 


0.267 


9.7 


10.0 


1918+26 


0.785 


4.1 


4.2 


1944+22 


1.334 


6.5 


6.7 


2027+37 


1.217 


4.8 


5.0 


2034+19 


2.075 


3.6 


3.7 


2035+36 


0.618 


3.6 


3.7 




WEISBERG 2004 - 0.430 G 


Hz 


Name B 


P 

[s] 


W50 (0.43 GHz) 

n 


W50 (1 GHz) 

n 


0045+33 


1.218 


5.6 


3.2 


1915+22 


0.426 


15.6 


14.3 


1922+20 


0.238 


17.0 


15.6 


1929+15 


0.314 


10.8 


9.9 


PARKES I - 1.374 GHz 


Name J 


P 

[s] 


W50 (1.374 GHz) 

n 


W50 (1 GHz) 

n 


0957-5432 


0.2035 


6.4 


6.6 


1049-5833 


2.2022 


5.4 


5.6 


1123-6102 


0.6402 


5.6 


5.8 


1144-6217 


0.8506 


4.2 


4.3 


1252-6314 


0.8233 


8.8 


9.1 


1348-6307 


0.9278 


11.0 


11.4 


1349-6130 


0.2594 


22.2 


22.9 
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PARKES IV - 1.374 GHz 


Name J 


P 

[s] 


14^50 (1-374 GHz) 

n 


1^50 (1 GHz) 

n 


0942-5552 


0.6644 


9. 


9.3 


1121-5444 


0.5358 


14.1 


14.6 


1157-6224 


0.4005 


11.7 


12.1 


1225-6408 


0.4196 


6.9 


7.1 


1326-5859 


0.4780 


5.7 


5.9 


1534-5408 


0.2897 


7.6 


7.8 


1703-4851 


1.3964 


3.4 


3.5 


1749-3002 


0.6099 


21. 


21.7 


1801-2920 


1.0818 


8. 


8.3 


1807-0847 


0.1637 


9. 


9.3 


1849-0636 


1.4514 


7.4 


7.6 


1857+0057 


0.5369 


4.9 


5.1 


19Q9-HQ0Q7 


1.0169 


2.8 


2.9 


19Q9-HQ254 


0.9899 


4.0 


4.1 


1910-0309 


0.5046 


5.4 


5.6 


PARKES VI - 1.374 GHz 


Name J 


P 

[s] 


Wm (1.374 GHz) 

n 


Wio (1 GHz) 

n 


1717-3953 


1.0855 


31. 


32.0 


1728-4028 


0.8663 


5.6 


5.8 


1822-0848 


0.8348 


10. 


10.3 


1827-0750 


0.2705 


12.5 


12.9 



